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/.  o  .  SUMMARY 

>  An  experimental  study, has  been- made  of  the  development  of  a  turbulent  boundary 
layer  after  an  interaction  with  a  normal  shock  wave  ( strong  enough  to  cause  a  local 
separation)  in  a  strong  adverse  pressure  gradient.  This  type  r'ftowffccurs  in  air  breathing 
engine  components  (e.g.  supersonic  intakes,  transonic  compressor  stages  and  supersonic 
diffusers),  is  poorly  understood  and  cannot  be  satisfactorily  predicted.  The  measurements, 
made  in  a  closed  duct,  extended  well  downstream  of  the  shock  wave  interaction.  Detailed 
results  for  the  flow  are  presented  and  used  to  support  two  major  conclusions,  firstly' jt  is 
shown  that  the  post  shock  adverse  pressure  gradient  has  a  large  effect  on  boundary  layer 
development  through  the  interaction  and  downstream  of  it.  Consequently  existing  results 
for  interactions  without  a  post  shock  pressure  gradient  should  not  be  used  as  a  model 
for  practical  flows  which  typically  have  strong  pressure  gradients  applied  downstream 
of  the  shock  wave.  The  second  conclusion  was  that  the  shock  wave  in  a  rectangular  duct 
produced  a  pronounced  stabilising  effect  on  the  downstream  flow.  Surface  flow  visualization 
suggests  that  this  stabilization  is  achieved  by  streamwise  vortices  shed  into  the  flow  from 
the  separated  region  formed  by  the  shock  wave.  The  implications  of  this  result  to  nominally 
two-dimensional  flow  situations  and  to  flows  with  weak  interactions  without  local 
separations,  is.  discussed. 
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NOTATION 


Symbol 

A 

CP 

c't 

G 

H 

M 

Mi 

Pr 

Ps 

pw 

Pr 

Pr0 

Re a 

Tv 

To 

Too 

Toi 

T„ 

T 

u 

Ui 

u 

UT 

X 

X 

y 

s 

s* 

So* 

Si* 


Definition 

constant  in  logarithmic  law  of  the  wall  (=  5  0) 
pressure  coefficient  referred  to  nozzle  exit  conditions 
skin  friction  coefficient  ( =  r0/^Ui2) 


form  factor  (S*/6) 

Mach  number 

local  free  stream  Mach  number 
molecular  Prandtl  number 
static  pressure 
wall  pressure 
total  pressure 

total  settling  chamber  pressure 

Reynolds  number  at  shock  wave  position 

wall  temperature 

local  stagnation  temperature 

settling  chamber  stagnation  temperature 

free  stream  stagnation  temperature 

stagnation  temperature  measured  by  probe  in  the  boundary  layer 

(rp-rw)/(r0i-rw) 


mean  velocity  in  x  direction 
free  stream  velocity 
ulVt 

velocity  scale  based  on  wall  shear  stress  ( =  (r0/p)*) 
distance  from  leading  edge  of  duct  floor 
distance  from  start  of  wall  pressure  rise 
vertical  distance  from  the  wall 
velocity  boundary  layer  total  thickness 


velocity  boundary  layer  displacement  thickness  ^  =  j*  ^1  — 

displacement  thickness  immediately  before  the  shock  wave  or  pressure  rise 
displacement  thickness  at  start  of  mixing  region 


normalised  enthalpy  thickness 


puTf 

PiUiToi 


A  1 


1.  INTRODUCTION 


This  report  presents  measurements  of  boundary  layer  development  in  a  strong  adverse 
pressure  gradient  after  an  interaction  with  a  normal  shock  wave  strong  enough  to  cause  a  local 
separation  in  the  boundary  layer.  This  type  of  flow  occurs  in  several  important  practical 
situations  such  as:  supersonic  intakes  (Stewart  &  Fisher  (1970)),  transonic  compressor  stages 
(Leblanc  &  Goethals  (1975)),  supersonic  diffusers  (Seddon  (1967))  and  also  on  transonic  aerofoils 
(Vidal,  Wittliff,  Catlin  &  Sheen  (1973),  Vidal  &  Catlin  (1977),  Inger  &  Mason  (1975),  Alstatt 
(1977),  Vidal  &  Kooi  (1976),  Alber,  Bacon,  Masson  &  Collins  (1973)).  It  is  very  poorly  predicted 
by  current  calculation  methods  (Vidal  &  Kooi  (1976)). 

It  would  have  been  satisfying  to  devise  an  experiment  that  was  relevant  to  all  these  flows. 

However,  the  present  experiment  achieved  this  aim  to  only  a  limited  extent.  Firstly  the  influence 
of  three  dimensional  effects  will  probably  be  quite  different  in  the  various  flows.  The  experimental 
layer  which  was  studied  developed  on  the  floor  of  a  rectangular  duct  and  therefore  the  results 
will  be  more  relevant  to  a  rectangular  supersonic  intake  than  for  example  a  transonic  aerofoil. 

Secondly,  the  experiment  involved  a  boundary  layer  developing  on  a  flat  plate  in  zero  pressure 
gradient  before  interacting  with  a  normal  shock  wave,  whereas  flow  on  transonic  aerofoils,  in 
supersonic  diffusers  and  transonic  compressors  have  a  favourable  pressure  gradient  before  the 
shock  wave,  and  some  supersonic  intakes  have  adverse  pressure  gradients  before  the  shock  wave. 

This  is  not,  however,  thought  to  be  an  important  difference,  because  of  the  overwhelming  changes 
to  a  boundary  layer  brought  about  by  impulsive  separation  at  a  shock  wave  and  the  vigorous 
mixing  which  precedes  reattachment.  Given  this  dominating  influence  of  the  shock  wave  on  the 
layer’s  behaviour,  the  initial  condition  of  the  layer  entering  the  interaction  will  have  a  negligible 
effect  on  its  downstream  development.  Differences  in  upstream  conditions  will  cause  layers  of 
different  thicknesses  to  be  generated  before  the  shock  wave,  but  this  just  means  an  increase  in 
the  size  but  not  the  structure  of  the  interaction  and  apparently  does  not,  affect  the  layer’s  down¬ 
stream  behaviour  (Little  (1967),  Vidal  &  Kooi  (1976)).  Therefore  the  absence  of  pressure  gradient 
before  the  shock  wave  is  argued  to  have  little  effect  on  the  generality  of  these  results. 

Previous  work  on  the  problem  has  been  of  three  types.  Firstly  there  have  been  engineering 
studies  with  the  aim  of  improving  the  performance  of  a  particular  piece  of  hardware  (e.g.  Brown, 

Nawrocki  &  Paley  (1968)).  Secondly  there  have  been  several  studies  of  normal  shock  wave- 
boundary  layer  interaction  with  no  (other)  pressure  gradients  applied  to  the  flow  (Seddon  (1967), 

Kooi  (1975),  Matear,  Bresh  &  Viegas  (1976),  East  (1976),  Sawyer,  East  &  Nash  (1977),  Abbiss, 

East,  Nash,  Parker,  Pike  &  Sawyer  (1976)).  Finally  in  two  cases  (Vidal  el  al.  (1973),  Padova, 

Falk  &  Wittliff  (1980))  flows  with  limited  post-shock  pressure  gradients  have  been  generated  in 
blowdown  Ludwieg  tubes  by  varying  the  porosity  of  a  perforated  nozzle.  However  the  pressure 
gradients  have  been  mild  and  their  effect  observed  over  a  very  short  distance  downstream  of  the 
shock  wave. 

There  exist  therefore,  experimental  data  on  the  effect  on  boundary  layer  development  of  a  ! 

normal  shock  wave  without  an  adverse  pressure  gradient  and  of  course  there  exist  numerous 
studies  of  subsonic  turbulent  boundary  layer  development  in  strong  adverse  pressure  gradients 
(see  Coles  &  Hirst  (1968)).  In  the  absence  of  data  on  the  combined  effect  of  the  two  agencies  it  is 
tempting  (see  Inger  (1975))  to  add  these  two  separate  effects  to  estimate  boundary  layer  develop¬ 
ment  in  an  adverse  pressure  gradient  after  an  interaction  with  a  normal  shock  wave.  However  t 

as  turbulent  boundary  layers  are  highly  non-linear  a  simple  addition  is  probably  incorrect. 

It  is  a  question  addressed  in  this  study  and  is  of  theoretical  as  well  as  practical  interest.  The 
theoretical  interest  arises  from  the  technique  of  studying  the  turbulent  mechanisms  within  a 
turbulent  boundary  layer  by  analysing  the  response  of  a  layer  to  discontinuous  or  step  changes  in 
boundary  conditions  (see  Clauser  (1956)).  A  normal  shock  wave  impinging  on  a  zero  pressure 
gradient  boundary  layer  can  be  considered  as  impulsively  applying  a  pressure  step  to  the  layer. 
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Seddon  (1967)  observed  the  layer’s  response  to  such  an  impulse  for  a  distance  of  fifty  (undisturbed) 
layer  thicknesses  which  is  further  than  any  other  worker  has  done.  His  results  showed  that  the 
effects  of  the  shock  wave  on  the  layer  were  still  evident  at  the  last  observation  station.  In  his 
flow  there  was  no  pressure  gradient  applied  after  the  shock  wave  and  it  may  be  that  the  response 
took  a  long  time  to  die  out  as  it  was  simply  decaying  in  the  absence  of  other  strong  forces.  It  is 
possible  in  the  flow  studied  here  that  the  strong  post  shock  pressure  gradient  will  dominate  the 
flow  and  quickly  erase  the  effects  of  the  shock  wave.  To  separate  the  effects  of  shock  wave  and 
pressure  gradient  we  require  reference  or  comparison  layers  that  have  one  of  the  agencies  (shock 
wave  or  pressure  gradient)  absent  but  otherwise  are  identical  to  the  main  flow  under  investigation. 
In  the  present  study,  therefore,  results  were  taken  in  a  reference  layer  that  developed  in  the  same 
adverse  pressure  gradient  as  the  main  flow  but  without  the  shock  wave.  Results  by  Seddon  (1967)1 
provide  the  other  reference  layer  with  a  normal  shock  wave  but  no  post  shock  adverse  pressure 
gradient.  With  these  three  sets  of  results  the  effect  on  boundary  layer  development  of  shock  wave, 
pressure  gradient  and  their  combination  can  be  evaluated. 


2.  EXPERIMENT 
2.1  Experimental  Rig 

The  experimental  rig  is  shown  diagrammatically  in  figure  1.  An  undried  air  supply  of 
4 -26  kg/s  was  heated  to  a  stagnation  temperature  of  403K  (to  prevent  subsequent  condensation 
shocks)  and  passed  into  a  settling  chamber  of  0-605  m  diameter  containing  six  wire  screens. 
The  flow  area  was  then  reduced  subsonically  in  a  duct  section  that  changed  from  circular  to 
rectangular.  The  area  reduction  ratio  from  the  settling  chamber  to  the  sonic  throat  was  40  to  1 
(details  are  given  in  Schofield  (1975)).  The  flow  was  then  accelerated  to  Mach  1-41*  giving 
an  exit  static  temperature  near  ambient.  Mach  1  -41  was  chosen  as  being  representative  of  flow 
velocities  before  shock  waves  in  supersonic  intakes,  and  because  previous  results  suggested  it 
would  generate  a  fairly  large  separation  bubble  in  the  turbulent  boundary  layer.  The  exit  size 
of  the  nozzle  was  99  mm  x  82  mm  high.  The  experimental  duct  entrance  was  also  rectangular 
but  smaller  in  both  dimensions  so  that  some  of  the  nozzle  flow  was  spilt  around  the  sides  of  the 
duct.  This  meant  that  none  of  the  nozzle  boundary  layers  entered  the  duct  and  that  the  heat  loss 
from  the  forward  portion  of  the  duct  (which  could  not  be  lagged)  was  minimized  by  having 
similar  flow  on  both  sides  of  each  wall.  The  duct  walls  downstream  of  the  nozzle  exit  were  well 
insulated  (see  figure  1)  to  try  to  produce  an  adiabatic  flow  inside  the  duct.  The  rig  was  made  from 
aluminium  castings  fitted  with  hardened  steel  leading  edges.  Continuous  rubber  seals  were  fitted 
between  the  castings  as  shown  in  figure  1. 

The  experimental  boundary  layer,  on  the  floor  of  the  duct,  started  from  a  leading  edge  in 
unsheared  flow,  had  a  natural  transition  to  a  turbulent  layer  and  developed  in  zero  pressure 
gradient  over  a  distance  of  0-24  m  before  interacting  with  the  normal  shock  wave.  Before  the 
interaction  the  layer  was  approximately  3-3  mm  thick  and  as  the  duct  width  was  85  mm,  there 
was  a  layer  thickness  to  duct  width  ratio  of  1/26. 

The  height  of  the  duct  entrance  was  47  mm  but  increased  with  distance  downstream  so 
that  the  post  shock  flow  was  subjected  to  an  adverse  pressure  gradient.  The  roof  was  made  of 
interlocking  plates  that  were  adjusted  to  generate  a  pressure  gradient  that  was  strong  but  did 
not  reseparate  the  layer  downstream  of  the  shock  wave.  The  sidewalls  of  the  duct  were  set  on 
slight  angles  to  the  centreline  of  the  floor  to  compensate  for  the  (calculated)  growth  in  the  sidewall 
boundary  layer  displacement  thickness.  Although  the  duct  width  thus  increased  (from  85  to 
100  mm)  the  total  layer  thickness  increased  more  rapidly  and  at  the  end  of  the  measurement 
region  the  layer  thickness  to  duct  width  ratio  was  approximately  1/3. 

Just  downstream  of  the  shock  wave,  fillets  were  attached  to  all  four  inner  corners  of  the  duct 
(see  figure  1)  to  reduce  secondary  flows.  The  thickness  of  the  fillets  increased  (from  zero)  with 

1  And  to  a  lesser  extent  by  Kooi  (1975),  East  (1976)  and  Abbiss  et  at.  (1976). 

*  Actually  1  -405±0-01  at  the  exit  plane,  see  Schofield  (1975). 
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distance  down  the  duct.  Surface  flow  visualization  indicated  that  the  duct  flow  (away  from  the 
shock  wave)  was  two  dimensional,  that  there  was  no  leakage  from  the  duct  and  that  the  boundary 
layers  on  the  other  three  walls  did  not  separate. 


2.2  Instrumentation 

The  duct  floor  contained  forty-six  small  diameter  static  pressure  tappings  spaced  along 
either  side  of  the  duct  centreline  (visible  in  figure  3 b).  In  addition  fifteen  20  mm  diameter 
instrumentation  holes  were  spaced  along  the  centreline  of  the  duct.  Into  these  holes  a  probe 
traversing  device  could  be  fitted  which  could  hold  a  variety  of  ‘gooseneck’  probes  (single  and 
double  total  head,  static  and  total  temperature).  When  not  in  use  these  holes  were  filled  with 
aluminium  plugs  set  flush  with  the  test  surface  and  sealed  with  ‘O-ring’  seals.  These  plugs  were 
instrumented  with  either  static  pressure  taps  or  copper-constantin  thermocouples  set  in  a  copper 
bead  flush  with  the  flow  surface. 

The  total  head  probes  used  near  the  wall  had  flattened  openings  measuring  (typically) 
1  mmxO- 12  mm  with  a  centre  of  pressure  estimated  to  be  at  0- 14  mm  from  the  lower  surface 
of  the  probe.  The  temperature  probe  was  a  standard  ventilated  stagnation  probe  with  a  copper- 
constantin  thermo-couple.  The  opening  of  this  probe  was  5  mm  x  0-38  mm  with  a  centre  of 
pressure  estimated  to  be  at  0  -35  mm  from  the  lower  surface  of  the  probe. 

All  thermocouples  were  referred  to  a  standard  ice  point  junction  and  the  voltages  generated 
were  read  with  a  digital  voltmeter  which  had  a  resolution  of  1  microvolt.  The  thermocouples, 
as  well  as  the  boundary  layer  and  the  settling  chamber  temperature  probe,  were  calibrated  against 
a  standard  mercury  and  glass  thermometer  in  water  over  the  range  0  to  100°C.  The  results  for 
each  thermocouple  were  fitted  with  a  third  order  polynomial  to  be  used  in  data  reduction.  As 
the  temperatures  measured  in  the  experiments  were  near  120°C  this  procedure  involved  the 
assumption  that  the  calibration  curves  could  be  accurately  extrapolated. 

The  static  pressure  tappings  were  connected,  using  small  diameter  tubing,  to  a  pressure 
scanning  valve  where  the  pressures  were  sensed  with  an  absolute  pressure  transducer.  Settling 
chamber  stagnation  pressure  and  atmospheric  pressure  were  read  in  the  same  way.  The  total 
pressure  probe  was  connected  to  one  side  of  a  differential  pressure  transducer  while  the  other 
side  was  connected  to  the  ganged  wall  static  pressure  tappings  either  side  of  the  probe  tip.  In 
this  way  the  dynamic  pressure  was  found  directly  with  one  transducer.  The  pressure  transducers 
were  periodically  calibrated  against  a  deadweight  tester  and  were  always  found  to  be  linear  to 
within  0  05%  of  full  scale.  The  sensitivity  for  a  given  transducer  changed  less  than  0  - 1  %  of  full 
scale  between  calibrations.  The  zero  pressure  reading  changed  much  more,  being  sensitive  to 
changes  in  temperature.  Consequently  the  zero  pressure  reading  was  taken  immediately  at  the 
conclusion  of  each  test  and  used  in  reducing  that  set  of  data. 

The  probe  traverse  incorporated  a  variable  speed  electric  motor  that  could  generate  steps 
as  small  as  0-008  mm.  The  probe  position  was  indicated  by  a  linear  distance  transducer  with  a 
travel  of  50  mm.  The  calibration  of  this  transducer  was  found  to  be  non-linear  and  affected  by 
temperature.  Consequently  it  was  calibrated  in  an  oven  over  the  temperature  range  it  attained 
on  the  rig.  The  two  hundred  calibration  points  were  fitted  using  a  cubic  spline  routine  that  was 
subsequently  employed  in  data  reduction.  The  temperature  of  the  transducer  on  the  rig  was 
measured  with  a  thermocouple. 

A  schlieren  system  allowed  continuous  viewing  of  the  shock  wave  region  through  the  windows 
in  the  duct  shown  in  figure  1.  The  system  also  allowed  spark  photographs  of  this  Schlieren 
pattern  to  be  taken. 


2.3  Uncertainty  Estimates 

The  uncertainty  of  the  experimentally  measured  quantities  were  estimated  as: 
y  =  ±0  008  mm 
T  =  ±0-8  K 
p  =  ±0-8%  of  reading 
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Estimates  of  uncertainty  for  derived  quantities  were  determined  using  the  method  of  Kline  & 
McClintock  (19S3)  and  are  given  below: 

u/Ui  =  ±1-2%;  ±2-5%  near  the  interaction  region;  ±4%  for  backflow  in  the  separation 
bubble. 

Pe/Pw  =  ±1  '7%  except  near  shockwaves 
T/Toi  =  ±0-3% 
c't  =  ±7% 


2.4  Procedure 

Firstly  the  ambient  pressure  was  measured  and  used  to  adjust  the  settling  chamber  blowing 
pressure  to  give  a  correctly  expanded  supersonic  jet  at  Mach  1  -41  from  the  nozzle.  The  rear  flap 
on  the  duct  was  then  slowly  closed  until  a  normal  shock  wave  appeared  on  the  Schlieren  screen. 
Flow  conditions  were  then  allowed  to  stabilize.  It  took  approximately  30  minutes  for  the  wall 
temperatures  at  the  rear  of  the  duct  to  become  constant.  The  shock  wave  position  was  then  adjusted 
to  coincide  with  a  marked  position  on  the  Schlieren  viewing  screen  where  it  was  maintained 
throughout  the  test.  The  wall  static  pressure  and  temperature  distributions  along  the  test  surface 
were  recorded  before  and  after  the  profile  data  were  taken.  At  each  point  in  a  pressure  profile 
the  probe  position,  probe  pressure  and  wall  pressure  were  recorded;  in  a  temperature  profile 
the  recordings  were:  probe  position,  probe  temperature  and  wall  temperature. 

After  measurements  in  the  main  flow  were  complete,  corresponding  measurements  in  the 
shockless  low  speed  reference  layer  were  made.  To  set  up  this  flow  the  settling  chamber  pressure 
was  reduced  so  that  the  nozzle  supplied  air  to  the  duct  at  Mach  0-61.  This  Mach  number  was 
lower  than  the  correct  post  shock  Mach  number  (0  ■  74)  behind  a  normal  shock  wave  at  Mach  1-41 
but  was  the  highest  that  could  be  generated  without  causing  some  supersonic  flow  with  shock 
waves  in  the  nozzle. 

On  the  basis  of  results  presented  below,  it  is  argued  that  this  difference  in  Mach  number 
is  insignificant:  it  will  be  seen  that  both  flows  have  very  similar  pressure  gradients  downstream 
of  the  normal  shock  position,  suggesting  no  significant  differences  due  to  compressibility. 

At  the  completion  of  these  measurements  the  back  flap  was  opened  fully  so  that  the  flow 
along  the  duct  accelerated  reaching  speeds  up  to  Mach  2-3  before  the  end  of  the  duct.  Again 
wall  pressures,  temperatures  and  profile  data  were  recorded.  These  supersonic  results  will  be 
reported  fully  in  a  separate  publication  but  some  of  the  results  are  used  in  this  report  to  check 
the  accuracy  of  the  temperature  profile  measurements. 


2.5  Data  Reduction 

Data  reduction  was  done  with  the  aid  of  a  computer  program  that  required  an  input  of  the 
raw  data  and  instrument  calibration  constants  for  each  test.  Profile  pressures1  were  converted 
to  Mach  numbers  using  the  isentropic  flow  relations.  If  a  temperature  profile  had  been  taken  for 
a  particular  station  then  it  was  interpolated  with  the  Mach  number  profile  to  give  a  velocity 
profile.  If  a  temperature  profile  was  not  available  at  that  station  then  the  simple  Crocco  (1932) 
relation, 

T=D  (1) 

was  used.  Comparisons  (given  later)  with  measured  temperature  profiles  show  that  equation(l) 
gave  satisfactory  correlations  for  zero  and  adverse  pressure  gradients.  In  favourable  pressure 
gradients  the  modified  Crocco  relationship 

T  =  0*  (2) 

was  found  to  be  more  representative  of  the  data. 

1  Using  the  probe  recovery  relationship  given  by  Winter  &  Gaudet  (1970). 
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The  data  reduction  program  calculated:  integral  boundary  layer  thicknesses,  static  tempera¬ 
ture  profiles  and  from  these,  viscosity  profiles  using  Sutherland’s  equation.  Density  profiles 
were  calculated  by  assuming  the  perfect  gas  law.  The  program  also  calculated  the  two 
dimensional  stream  function 


and  a  normalised  enthalpy  thickness 
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Reduced  data  is  tabulated  in  Table  1.  Layer  S  denotes  the  main  flow  with  the  shock  wave 
boundary  layer  interaction  followed  by  the  adverse  pressure  gradient.  Layer  L  denotes  the 
reference  flow  without  a  shock  wave  but  with  the  adverse  pressure  gradient. 

Skin  friction  coefficients  were  determined  using  Clauser’s  method  (Clauser  (1954))  which  is 
based  on  writing  the  logarithmic  law  of  the  wall  in  the  form. 
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and  comparing  data  points  on  u/Ui,  yUJv  axes  with  a  family  of  straight  lines,  the  slope  and 
position  of  which  depends  on  uTIU\(=\Jct'l2).  (The  profiles  are  represented  in  the  Appendix.) 
As  it  was  calculated  in  the  data  reduction  process,  the  local  kinematic  viscosity  was  used  in 
calculating  yUi/v  values.  Otherwise  no  allowance  was  made  for  compressibility.  Only  the  two 
profiles  upstream  of  the  shockwave  in  series  S  would  be  significantly  affected  by  compressibility. 
At  these  Mach  numbers  the  differences  in  estimated  skin  friction  between  an  incompressible 
analysis  and  a  compressible  one  is  of  order  10%.  As  this  variation  is  about  the  difference  in  values 
of  skin  friction  given  by  different  compressible  methods  and  as  the  actual  values  for  these  profiles 
were  unimportant  in  the  analysis  of  this  flow,  only  the  incompressible  values  have  been  quoted. 


3.  RESULTS 
3.1  Flow  Visualization 

Figure  2  shows  a  Schlieren  spark  photograph  of  the  shock  wave.  It  is  similar  to  a  photograph 
presented  by  Seddon  (1967)  for  his  interaction  with  a  zero  pressure  gradient  boundary  layer. 
The  normal  shock  wave  bifurcates  about  6  to  8  layer  thicknesses  above  the  wall,  into  a  pair  of 
forward  and  rearward  facing  oblique  shock  waves.  A  fine  line  springing  from  the  bifurcation 
point  indicates  the  existence  of  a  vortex  sheet.  Such  a  shock  structure  suggests  that  a  large 
separated  region  exists  under  the  shock  wave.  The  surface  flow  pattern  (figures  3a,  3b)  confirms 
that  there  is  a  large  area  of  separated  flow  and  shows  that  it  is  strongly  three  dimensional. 
Figure  3c  presents  an  interpretation  of  the  wall  flow  patterns.  It  shows  the  flow  separating,  under 
the  shock  wave,  along  lines  (separatrices)  emanating  from  a  separation  saddle  point  and  curving 
slightly  downstream  as  they  approach  the  duct  sidewalls.  Directly  downstream  of  the  separation 
saddle  is  a  reattachment  saddle  into  which  flow  two  reattachment  lines.  On  each  side  of  a  line 
joining  the  two  saddles  the  separated  flow  winds  up  into  two  counter-rotating  focii. 

A  very  similar  pattern  (and  interpretation)  is  given  in  Green  (1969)  for  an  interaction  of  a 
turbulent  boundary  layer  with  a  strong  oblique  shock  wave  in  a  square  duct.  Very  similar 
separated  regions  are  formed  behind  blocks  attached  to  the  surface  in  a  boundary  layer  as  shown 
in  figure  4  (see  also  Castro  &  Robins  (1977)).  Hunt,  Abell,  Peterka  &  Woo  (1978)  made  an 
experimental  and  analytical  study  of  these  separated  regions.  Two  of  their  conclusions  that  are 
important  to  this  work  are:  there  are  no  closed  surfaces  around  the  separated  flows  and  that 
counter-rotating  vortices  are  shed  into  the  downstream  flow  fi .  >  the  separated  region.  Supporting 
evidence  for  this  last  conclusion  is  given  by  Colmer  (1970)  who  measured,  at  distances  well 
downstream  of  a  rectangular  aircraft  hanger,  significant  upfiows  on  the  centreline  of  flow, 
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suggesting  the  presence  of  streamwise  vortices.  Other  indirect  evidence  is  presented  by  Schofield  & 
Logan  (1983)  who  found  for  the  two  cases  shown  in  figure  4  that  the  mean  centreline  flow 
recovered  more  quickly  behind  the  narrower  of  the  two  blocks.  The  explanation  offered  is  that 

the  narrow  block  the  shed  vortices  are  closer  to  the  centreline  and  thus  their  mixing  action 
is  more  effective  on  the  centreline  flow  giving  a  quicker  recovery. 

It  is  theoretically  possible  that  in  the  present  flow  vortices  shed  from  the  two  focii  join 
up  and  are  contained  within  a  closed  separation  bubble.  In  view  of  the  preceding  evidence  it 
seems,  however,  unlikely.  The  existence  of  longitudinal  vortices  in  the  present  flow  offers  the 
most  plausible  explanation  of  the  results  presented  in  this  report. 

The  theory  of  flow  around  critical  points  (saddles,  focii  and  nodes)  in  viscous  flows  has  been 
developed  considerably  in  the  last  few  years  (Perry  &  Fairlie  (1974),  Hunt  et  al.  (1978),  Perry, 
Chong  &  Lim  (1982)).  One  aspect  of  the  work  presented  in  these  papers  is  the  supply  of  mass  to 
vortices  as  they  roll  up  and  flow  downstream.  In  the  present  case,  mass  would  have  to  flow 
into  the  separated  region  to  supply  the  vortices  shed  into  the  downstream  flow.  The  separated 
region  cannot  therefore  be  closed,  but  must  have  stream  tubes  correcting  it  with  the  upstream 
flow.  Thus,  because  open  streamtubes  may  exist  dose  to  the  wall,  the  separation  and  reattach¬ 
ment  lines  have  not  been  joined  in  figure  3c  although  the  oil  flow  pattern  (which  is  fairly  indistinct 
in  this  region)  could  be  interpreted  as  supporting  their  junction. 

3.2  Pressure  Field 

The  wall  pressure  distributions  measured  in  the  two  present  flows  are  shown  in  figure  5  and 
the  corresponding  free  stream  Mach  number  distributions  in  figure  6.  Although  the  Mach  number 
after  the  shock  wave  in  layer  5  was  higher  than  the  nozzle  Mach  number  in  layer  L,  differences 
in  turbulence  structure  due  to  compressibility  would  have  been  very  small.  Importantly,  pressure 
gradients  of  the  two  flows  appear  very  similar  (see  figures  7  and  8).  Since  pressure  gradient  is 
the  dominating  factor  in  boundary  layer  development  it  is  argued  that  layer  L  is  a  good  reference 
layer. 

In  both  separating  flows  investigated  by  Simpson  et  al.  (1977  and  1981)  separation  occurred 
just  before  the  pressure  gradient  dropped  to  a  near  zero  constant  value.  In  the  present  results 
separation  occurs  before  a  constant  value  of  pressure  gradient  is  attained.  Perry  &  Fairlie  (1974) 
have  presented  analysis  and  experimental  results  which  suggest  that  separation  and  reattachment 
points  coincide  with  the  pressure  gradient  becoming  (locally)  zero.  As  shown  in  figure  8  the 
present  results  do  not  support  this  proposal. 

The  shape  of  the  initial  wall  pressure  rise  under  the  shock  wave  appears  to  be  governed 
mainly  by  the  shock  strength.  Figure  9  shows  that  the  wall  pressure  rise  in  the  early  stages  of  the 
interaction  in  the  present  flow  is  similar  to  that  in  Seddon’s  flow  which  did  not  have  a  post-shock 
pressure  gradient.  To  compare  the  results  the  downstream  distance  has  been  non-dimensionalized 
with  the  undisturbed  boundary  layer  displacement  thickness.  It  is  seen  later  that  in  the  present 
flow,  the  post-shock  adverse  pressure  gradient  produces  large  changes  to  the  flow  structure 
but  apparently  it  has  little  effect  on  the  shape  of  the  initial  pressure  rise.  This  is  consistent  with 
the  concept  of  a  ‘free  interaction’  proposed  by  Chapman,  Kuehn  &  Larsen  (1957)  in  which, 
provided  a  shock  wave  induced  separation,  the  flow  upstream  of  the  separation  point  was  not 
muw  affected  by  downstream  conditions. 

The  Prandtl  boundary  layer  approximation  leads  to  0  which  is  a  very  good  approxi¬ 

mation  for  layer  L  (the  flow  without  a  shock  wave)  as  shown  in  figure  10.  Near  a  shock  wave 
however,  it  is  certain  that  dpa/dy  cannot  be  zero  as  the  streamwise  pressure  distribution  changes 
from  the  wall  distribution  shown  in  figure  5  to  a  step  distribution  in  the  flow  above  the  shock 
bifurcation  point.  The  surface  flow  pattern  of  figure  3  suggests  that  the  static  pressure  will  also 
vary  laterally  across  the  flow.  The  static  pressure  field  near  the  separation  bubble  varies,  therefore, 
in  all  three  directions. 

In  the  present  work  the  static  pressure  was  assumed  to  be  constant  through  the  boundary 
layer.  A  few  static  pressure  profiles  were  recorded  (figure  1 1)  and  they  show  this  to  be  an  acceptable 
approximation  for  flow  near  reattachment  and  downstream  from  it.  Upstream  of  reattachment, 
near  separation  and  the  shock  waves,  is  a  most  difficult  area  to  obtain  accurate  static  pressure 
measurements  in  the  flow  and  it  was  not  attempted  here.  However  the  magnitude  or  even  the 
sign  of  departure  from  the  constant  pressure  assumption  is  far  from  clear.  Seddon’s  limited 
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measurements  in  this  region  show  no  static  pressure  variation  within  the  boundary  layer  but 
significant  variation  outside  it  up  to  the  height  of  the  bifurcation.  Kooi  (1975)  measured  four 
static  profiles  near  the  shock  wave  but  all  were  quite  different  from  Seddon’s  result  as  they  showed 
significant  variation  within  the  boundary  layer.  To  further  confuse  the  issue  Kooi's  profiles 
were  radically  different  in  shape  between  themselves,  even  those  spaced  only  0-88  apart.  This 
conflicting  data  suggests  that  the  assumption  of  constant  static  pressure  is  as  good  as  any  at 
present.  This  uncertainty  is  reflected  in  the  larger  uncertainty  estimate  for  ujUi  measurements 
in  the  interaction  region  (Section  2.3). 


3.3  Velocity  Field 

Velocity  profiles  were  measured  only  on  the  centreline  of  the  flow.  In  the  interaction  region 
profiles  off  the  centreline  will  vary  significantly.  Downstream  of  the  interaction,  the  surface 
flow  tests  suggest  that  the  centreline  should  be  representative  of  the  flow  over  most  of  the  width. 
The  velocity  and  density1  profiles  have  been  used  to  generate  the  function 

<l>  =  J*V<a!y. 

For  two  dimensional  flow  i/>  is  the  stream  function.  Areas  of  the  interaction  flow  are,  however 
highly  three  dimensional  and  in  these  regions  4>  will  not  give  streamlines.  However  the  present 
mean  velocity  measurements  were  taken  on  the  duct  centreline  and  here  flow  symmetry  would 
require  three  dimensionality  or  crossflows  to  be  negligible.  On  the  centreline  then  > p  will  probably 
approximate  the  stream  function.  Figure  12  shows  the  centreline  streamlines  with  the  (Mach) 
profiles  used  to  calculate  them.  An  average  shock  system  taken  from  the  spark  schlieren  photo¬ 
graphs  is  also  shown. 

The  data  in  this  diagram  is  internally  consistent,  the  streamline  flow  angles  show  good  agree¬ 
ment  with  the  flow  deflections  of  the  shock  system  and  the  values  of  </>  given  by  different  profiles 
form  a  consistent  streamline  pattern.2  However  there  are  four  factors  that  could  affect  the  accuracy 
of  this  diagram : 

(i)  static  pressure  variations  through  the  boundary  layer, 

(ii)  As  figure  36  shows,  the  line  of  symmetry  of  the  separated  flow  does  not  quite  co¬ 
incide  with  the  duct  centreline  and  thus  crossflows  on  the  duct  centreline  may  not 
be  negligible. 

(iii)  velocity  vectors  have  been  assumed  to  be  normal  to  the  pitot  probe  opening  whereas 
Abbiss  et  al.  (1976)  showed  that  near  the  shock  wave,  flow  angles  could  be  as  large 
as  8°  to  the  x  axis, 

(iv)  the  reversed  profiles  were  measured  with  (small)  pitot  probes  and  Simpson  (1976) 
claims  that  they  cannot  accurately  measure  backflow  with  its  intermittent  flow 
reversals  and  consequent  large  changes  in  velocity  magnitude  and  direction. 

The  effects  of  (ii)  and  (iii)  can  be  estimated  and  are  small.  The  magnitude  of  (i)  cannot  be 
confidently  assessed  but  is  unlikely  to  be  significant  over  the  entire  flow  map.  Finally  the  magni¬ 
tudes  of  the  reversed  flow  inferred  from  the  pitot  tubes  was  similar  to  those  found  by  Simpson 
et  al.  (1977,  1981)  using  laser  anemometry.  Figure  12  therefore,  is  probably  a  fairly  accurate 
drawing  of  the  centreline  streamline  pattern. 

1  Density  was  calculated  by  the  data  reduction  program  from  the  perfect  gas  law.  using  either 
the  calculated  or  measured  static  temperature  profile  and  assuming  constant  static  pressure. 

2  The  separation  point  has  been  shown  connected  to  the  reattachment  point  although  it  is  quite 
possible  that  an  ‘alleyway’  exists  in  the  vertical  plane  to  supply  mass  into  the  separated  region. 
However,  the  streamline  plots  in  the  figure  are  not  accurate  enough  to  determine  such  detail 
(see  main  text). 
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Figure  13  compares  a  non-dimensionalized  version  of  this  streamline  pattern  with  the  one 
by  Seddon  for  a  flow  with  no  post  shock  pressure  gradient.  There  are  substantial  differences 
between  them  which  can  be  explained  by  the  effect  the  post  shock  adverse  pressure  gradient  has 
in  delaying  reattachment  and  hence  creating  a  larger  separation  region  than  in  Seddon’s  flow. 
In  the  present  flow  the  rapid  increase  in  displacement  thickness  caused  by  the  larger  separated 
region  requires  larger  streamline  deflections.  Consequently  the  shock  wave  system  is  tipped 
forwards  and  upwards  away  from  the  wall  compared  with  the  pattern  measured  by  Seddon. 
Mach  number  profiles  suggest  that  a  large  supersonic  tongue  is  formed  behind  the  shock  wave. 
Seddon’s  measurements  also  suggested  a  supersonic  tongue  but  his  flow  had  a  higher  incident 
Mach  number  (1-47).  Flows  at  the  same  Mach  number  (1-4)  without  a  post  shock  pressure 
gradient  (East  (1976),  Kooi  (1975))  did  not  apparently  contain  supersonic  tongues.  Apparent 
inconsistencies  between  these  results  are  impossible  to  resolve  with  the  current  evidence.  All 
results  involve  some  uncertainty  in  the  Mach  number  profiles  near  the  shock  wave,  mainly  due 
to  inaccuracies  in  measuring  the  static  pressure.  Within  the  supersonic  tongues  the  Mach  numbers 
exceed  unity  by  small  amounts  and  small  differences  in  static  pressure  are  sufficient  for  a  super¬ 
sonic  tongue  to  appear  or  disappear  in  the  results. 


3.3.1  Displacement  Thickness 

The  distribution  of  displacement  thickness  in  the  interaction  region  is  compared  with  those 
for  flows  without  post  shock  pressure  gradients  in  figure  14.  The  three  comparison  flows  have 
similar  Mach  and  Reynolds  numbers  but  have  much  smaller  growths  in  displacement  thickness 
through  the  interaction  than  the  present  layer.  After  reattachment  the  displacement  thickness 
of  the  layers  without  a  post  shock  pressure  gradient  continues  to  decrease.  This  is  because  the 
layer  is  relaxing  back  to  zero  pressure  gradient  flow  conditions  and  velocities  near  the  wall  increase 
rapidly,  causing  the  displacement  thickness  to  decrease.  In  the  present  flow  however,  the  layer 
after  reattachment  continues  to  be  subjected  to  a  severe  adverse  pressure  gradient  which  causes 
high  rates  of  mass  entrainment  into  the  layer.  Although  flow  acceleration  near  the  wall  again 
would  tend  to  decrease  the  displacement  thickness,  the  entrainment  rate  of  the  layer  is  so  large 
that  in  this  case  the  displacement  thickness  increases  quite  rapidly  after  reattachment. 

Figure  15  compares  the  early  displacement  thickness  growth  for  the  present  layers  S  and  L. 
The  comparison  shows  that  in  the  region  where  the  pressure  gradients  for  the  two  layers  are 
similar  (i.e.  after  the  reattachment  position)  the  displacement  thickness  growth  rates  for  them 
are  very  similar.  The  growth  rates  are  high  while  the  local  pressure  gradient  is  high  but  when 
the  pressure  gradients  drop  (see  figure  8)  the  layers  go  through  identical  relaxations  where  the 
velocity  near  the  wall  increases  and  the  displacement  thickness  therefore  decreases.  At  the 
completion  of  this  process  the  displacement  thicknesses  begin  to  increase  again  by  entrainment. 
The  rates  are  closely  similar,  which  might  be  expected  as  the  two  layers  are  developing  in  similar 
pressure  gradients.  This  figure  would  suggest  therefore  that  the  effect  of  the  shock  wave  on  dis¬ 
placement  thickness  is  limited  to  a  simple  constant  addition  to  the  thickness.  It  might  appear 
therefore  that  the  effects  of  shock  wave  and  pressure  gradient  are  independent,  with  the  shock 
wave  having  only  a  local  effect  on  the  layer  (shown  in  figure  1 5)  and  the  pressure  gradient  deter¬ 
mining  boundary  layer  development  downstream  of  the  shock  wave.  However  this  is  not  the  case. 
Figure  16  plots  the  displacement  thickness  growth  for  the  two  layers  through  the  whole  down¬ 
stream  mixing  region.  To  give  a  fair  comparison  of  growth  rates  in  this  mixing  region,  the 
displacement  thickness  has  been  non-dimensionalized  with  its  value  at  the  start  of  the  region 
(x  =  0-403  m);  in  this  way  the  difference  in  displacement  thickness  due  to  the  shock  wave  in 
layer  S  is  removed  from  the  downstream  comparison.  The  figure  shows  that  the  displacement 
thickness  grows  much  faster  for  the  layer  without  the  shock  wave  and  would  probably  separate 
before  the  layer  which  has  interacted  with  a  shock  wave. 

The  explanation  proposed  here  for  this  difference  between  the  two  layers  is  that  vortices 
are  shed  from  the  separated  region  under  the  shock  wave  and  have  an  important  effect  on  flow 
downstream  of  the  high  pressure  gradient  region  by  increasing  mixing  in  the  layer,  stabilizing 
it  and  slowing  its  growth.  This  conclusion  is  supported  by  analysis  of  other  boundary  layer 
parameters  in  the  following  sections. 
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3.3.2  Boundary  Layer  Shape  Factors 

In  figure  17,  the  standard  boundary  layer  shape  factor  (8 */#)  is  compared  with  Seddon’s 
results  for  no  post  shock  pressure  gradient.  The  main  difference  between  the  results  is  in  the 
separation  region  and  is  due  to  the  enlargement  and  extension  of  the  separated  flow  by  the  adverse 
pressure  gradient  in  layer  S.  After  reattachment  the  results  are  similar  in  this  interaction  region. 

Figure  18  compares  the  shape  factor  variation  with  that  of  the  shockless  reference  layer 
through  the  downstream  mixing  region.  These  results  again  indicate  that  the  shockless  layer 
will  separate  first  as  it  has  a  larger  value  of  H  throughout  the  mixing  region  and  is  increasing 
quite  rapidly  at  the  last  observation. 

Similar  conclusions  can  be  drawn  from  the  variations  of  two  other  shape  factors:  Clauser’s 
shape  factor,1  (G,  figure  19)  and  Cole's  wake  strength  parameter2  (FI,  figure  20). 


3.3.3  Skin  Friction 

Comparisons  of  skin  friction  coefficient  for  the  two  present  layers  (figure  21)  show  that 
downstream  of  the  interaction  the  shockless  layer  is  moving  towards  separation  while  the  results 
of  the  layer  with  a  shockwave  are  nearly  constant  with  distance. 

In  the  interaction  region  where  shockwave  forces  dominate,  comparisons  with  flows  that 
do  not  have  post  shock  pressure  gradients  are  interesting.  Some  results  given  by  Seddon  (figure  22) 
show  that  the  skin  friction  coefficient  has  overshot  the  undisturbed  value  while  other  results  by 
Kooi  (1975)  suggest  that  they  may  also  overshoot  the  undisturbed  value  further  downstream. 
Figure  22  shows  that  the  present  results  don’t  do  this.  A  simple  explanation  for  this  difference 
is  that  in  all  the  flows  longitudinal  vortices  are  generated  at  the  interaction  and  have  a  strong 
stabilizing  effect  on  the  downstream  flow.  In  the  absence  of  other  forces  these  vortices  quickly 
increase  the  skin  friction  coefficient  of  the  flows  above  their  upstream  undisturbed  values.  In 
the  present  flow  however  the  effect  of  the  vortices  are  counteracted  in  part  by  the  destabilizing 
effect  of  the  post  shock  adverse  pressure  gradient.  This  interaction  results  in  the  skin  friction 
stabilizing  at  a  lower  value  than  the  upstream  undisturbed  value. 

It  seems  from  figure  21  that  the  vortex  mixing  persists  for  a  long  distance  downstream. 
This  long  persistence  of  the  mixing  vortices  is  important  in  considering  the  question  raised  by 
Seddon:  what  length  of  flow  development  is  required  for  the  layer  to  return  to  its  undisturbed 
state?  The  above  discussion  implies  that  for  layers  with  or  without  post  shock  pressure  gradients, 
the  mixing  vortices  initiated  at  the  separation  bubble  would  have  to  decay  before  the  relaxation 
process  can  commence.  The  present  results  show  no  indication  of  such  a  decay  at  the  last 
observation  (which  corresponds  to  an  x/S*  value  of  4000)  and  it  seems  likely  therefore  that  a 
return  to  an  ‘undisturbed  state'  will  require  a  distance  so  large  that  it  would  be  of  little  practical 
interest. 


3.4  Temperature  Field 
3.4.1  Wall  Temperatures 

The  centreline  wall  temperature  variation  for  the  shock  wave  layer  (S)  is  compared  with 
the  theoretical  adiabatic  wall  temperature  distribution3  for  zero  pressure  gradient  flows  in  figure  23. 
As  the  measured  temperature  distribution  is  lower  in  value  than  the  adiabatic  distribution  there 
is  a  net  heat  flow  from  the  plate  to  its  lower  temperature  surroundings.  Near  the  reattachment 
point  the  wall  temperature  approaches  the  stagnation  temperature  but  this  is  followed  down¬ 
stream  by  a  local  minimum  in  wall  temperature  which  coincides  with  the  region  of  large 
accelerations  in  the  flow  near  the  wall.  Downstream  of  this  region  the  wall  temperature  is  main¬ 
tained  at  a  little  less  than  the  (zero  pressure  gradient)  adiabatic  distribution. 

1  Clauser  (1954). 

2  Coles  (1956). 

3  Using  a  wall  recovery  factor  of  Prl/3. 
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Figure  24  shows  the  normalized  enthalpy  thickness  as  a  proportion  of  total  layer  thickness 
along  the  flow.  The  normalized  enthalpy  thickness  of  a  boundary  layer  is  that  thickness  of  flow 
at  free  stream  conditions  that  would  account  for  the  difference  between  local  and  free  stream 
enthalpy  integrated  across  the  boundary  layer.  It  is  a  (total)  heat  energy  displacement  thickness 
and  is  positive  for  heat  flowing  into  the  layer  and  negative  for  heat  flowing  out  of  the  layer. 
Figure  23  shows  that  it  is  mainly  negative  along  this  flow  but  in  all  cases  is  small  amounting  to 
less  than  0-5%  of  the  total  heat  in  the  boundary  layer. 

The  heat  loss  from  the  plate  in  layer  L  is  greater  as  indicated  by  The  <"  r. 

distribution1  shown  in  figure  25.  In  this  case  the  temperature  distribution  varies  slowly  and 
continuously  with  distance  and  is  much  further  below  the  zero  pressure  gradient  wall  adiabatic 
temperature  distribution.  This  is  to  be  expected  as  the  temperature  of  the  air  supplied  to  the 
experimental  duct  was  the  same  for  both  layer  L  and  S  but  the  mass  flow  and  hence  heat  supplied 
for  layer  L  was  only  38%  of  that  supplied  for  layer  S'.  The  duct  heat  loss'-  for  layer  L  was  thus  a 
!  r  oroportion  of  the  heat  supplied  to  the  duct  and  resulted  in  the  larger  temperature  drops 
shown  in  figure  25. 


3.4.2  Temperature  Profiles 

The  temperature  profiles  measured  here  in  strong  adverse  pressure  gradients  are  unusual, 
particularly  those  through  or  near  the  separated  region.  As  the  temperature  differences  between 
the  flow  and  the  wall  were  small  and  the  probe  was  a  simple  shielded  thermocouple  it  was  thought 
desirable  to  check  the  reliability  of  the  instrumentation  and  measuring  techniques.  This  was 
done  by  measuring  temperature  profiles  in  a  shockfree  adiabatic  favourable  pressure  gradient 
flow  and  comparing  them  with  reliable  data  that  exists  for  these  flow  conditions.  The  flow  was 
generated  by  opening  the  rear  flap  of  the  duct  to  give  an  accelerating  supersonic  flow  with  Mach 
numbers  up  to  2  -3.  The  wall  temperature  results,  figure  26,  shows  that  the  centreline  distribution 
was  ciom.  .v.  v_..  _ —  "'aHienl)  adiahr.'ic  conditions  The  temperature  profiles  (figure  27) 

may  be  compared  with  those  of  Meier  &  tun  ia  (!/<•,  ' c * ' 

in  a  supersonic  nozzle  using  a  most  sophisticated  probe  that  drew  air  through  a  range  oi  um... 
orifices.  Both  the  present  flow  and  that  of  Meier  &  Rotta  were  closely  adiabatic  and  had  mild 
favourable  pressure  gradients.  The  data  are  thus  directly  comparable. 

Both  sets  of  temperature  profiles  (the  present  ones  and  those  of  Meier  &  Rotta)  have  a 
sinuous  shape  that  crosses  the  Crocco  line  (T  =  0)  before  overshooting  the  free  stream  stag¬ 
nation  temperature  (7"  =  1).  The  size  of  the  temperature  overshoots  were  slightly  different  in 
the  two  layers.  For  adiabatic  flow,  conservation  of  energy  requires  that  the  average  stagnation 
;„-rr)sS  the  layer  should  equal  the  free  stream  enthalpy.  This  is  obviously  not  fulfilled 
in  the  present  prouics  uS  the  temperature  o are  too  small  to  compensate  for  the 

temperature  deficits  near  the  wall.  However  they  are  consistent  . .  K--»t  flow-  out  of 

the  layer.  As  the  Meier  &  Rotta  flow  was  closer  to  adiabatic  conditions  than  tne 

their  profiles  had  larger  overshoots.  Overall  then  the  consistency  of  these  results  gave  some 

confidence  in  the  present  temperature  measurements. 

The  temperature  profiles  in  the  shock  free  adverse  pressure  gradient  flow  (layer  L)  are  quite 
different  from  the  favourable  pressure  gradient  results.  Figure  28  shows  that  while  these  profiles 
still  have  a  sinuous  shape  they  do  not  cross  the  Crocco  line  ( T  =  (?)  but  lie  between  it  and  the 
free  stream  stagnation  line.  None  of  the  profiles  show  an  overshoot  which  is  consistent  with  a 

rively  large  heat  flow  out  of  this  layer. 

As  rnigh*  be  exnected.  the  temperature  profiles  for  the  shockwave  layer  ( S )  differ  from  the 
previous  results  only  m  the  vicinity  <*•  m*»  i»ier.«non.  i-icu.c  2-.’  «>•  »f  flu. 

interaction  region  (x>0-4  m),  the  profiles  are  similar  to  those  in  Uyu  L 

they  have  overshoots  in  total  temperature  as  expected  for  near  adiabatic  conditions. 

1  And  the  values  of  t/8,  see  table  1. 

2  As  the  density  is  lower  in  layer  L,  the  heat  transfer  coefficient  from  the  flow  to  the  duct  wall 
will  be  lower  and  thus  the  overall  heat  transfer  from  the  duct  will  be  slightly  less.  This  however 
will  be  a  small  effect  compared  with  the  large  difference  in  duct  mass  flow. 
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The  profiles  around  the  separated  region  are  quite  unusual.  As  yet  there  is  no  theory  that 
can  suggest  what  mechanisms  in  the  boundary  layer  cause  such  distributions  in  total  temperature. 
The  best  that  can  be  done  is  to  try  and  correlate  the  flow  behaviour  at  each  station  with  the  total 
temperature  distributions.  The  profile  at  x  =  0-249  m  passes  through  the  recirculating  fluid1 
and  displays  a  nearly  constant  total  temperature  across  the  layer.  The  outer  part  of  the  profile 
suggests  that  the  upstream  wall  fluid  lifted  over  the  backflow  has  been  fully  mixed.  This  is 
supported  by  the  results  of  Simpson  et  al.  (1981)  who  made  Reynolds  stress  measurements  after 
detachment  which  implied  very  high  rates  of  mixing  just  outside  the  reversed  flow  region.  The 
profile  at  x  =  0-276m  has  a  large  overshoot  and  is  nearly  adiabatic  (see  figure  24).  It  was 
positioned  very  close  to  the  reattachment  saddle  where  the  flow  is  strongly  decelerating  and  is 
similar  to  a  stagnation  flow  with  a  wall  temperature  that  approaches  the  adiabatic  value.  In 
contrast  the  profile  at  x  =  0-314  m  has  a  low  average  total  temperature.  It  was  positioned  in  a 
region  where  the  wall  flow  was  undergoing  rapid  acceleration. 


4.  DISCUSSION  AND  CONCLUSION 

The  major  conclusion  of  this  experiment  is  that  the  effects  of  the  normal  shock  wave  and  the 
adverse  pressure  gradient  did  not  combine  to  produce  a  boundary  layer  near  separation  but 
rather  that  the  interaction  acted  as  a  stabilizing  agency  for  the  boundary  layer  in  the  adverse 
pressure  gradient.  This  result  will  be  relevant  to  flows  in  most  rectangular  supersonic  air  intakes 
and  diffusers.  Its  relevance  to  axisymmetric  air  intakes,  transonic  aerofoils  and  transonic  com¬ 
pressor  stages  is  more  doubtful.  While  it  is  theoretically  possible  that  these  flows  are  two 
dimensional  with  closed  separation  bubbles,  in  practice  however  the  flows  will  be  three 
dimensional  to  varying  degrees  and  may  even  develop  cellular  flow  patterns  shedding  vorticity 
downstream.  Although  the  level  of  this  vorticity  is  likely  to  be  considerably  less  than  that  in  the 
present  experiment,  it  may  have  a  significant  effect  on  the  downstream  layer.  However,  there  is 
another  possible  source  of  vorticity  which  could  be  important  in  these  nominally  two  dimensional 
flows.  It  is  likely  that  Taylor -Goertler  vortices  would  be  generated  in  the  flow  over  the  rear  of 
the  separation  bubble  as  a  consequence  of  the  change  in  angular  velocity  in  this  region.  These 
speculations  are  somewhat  supported  by  the  work  of  Livesey  &  Odukwe  (1974),  Odukwe  & 
Livesey  (1974)  and  Kamal  &  Livesey  <  1977).  These  workers  have  measured  the  performance  of 
a  range  of  subsonic  axisymmetric  diffusers  preceded  by  a  long  inlet  pipe.  They  found  that  the 
diffusers  gave  higher  overall  pressure  recovery  and  lower  outlet  distortion  (two  factors  closely 
related  to  boundary  layer  separation)  if  the  inlet  pipe  flows  contained  normal  shock  waves. 
The  origin  of  this  stabilization  was  attributed  to  the  (measured)  increased  turbulence  level. 
The  differences  in  diffuser  performance  were  large,  suggesting  an  intense  distributed  mixing 
agency  such  as  longitudinal  vortices. 

The  present  results  are  probably  not  relevant  to  the  case  where  the  shock  wave  is  too  weak 
to  separate  the  boundary  layer.  Here  the  adverse  effects  of  shock  wave  ,.nd  pressure  gradient 
on  boundary  layer  development  may  be  additive  in  some  way.  Paradoxically,  component 
performance  will  be  enhanced  in  some  cases  by  the  presence  of  strong  normal  shock  wave  with 
its  local  separation,  as  this  has  the  same  effect  as  a  pair  of  vortex  generators  at  the  start  of  the 
subsonic  pressure  rise. 

Another  conclusion  of  this  work  is  that  the  post  shock  pressure  gradient,  as  well  as  having 
a  large  effect  on  the  downstream  boundary  layer  development,  has  a  gross  effect  on  the  structure 
of  the  interaction  region  itself.  Compared  with  the  zero  pressure  gradient  case,  the  adverse 
pressure  gradient  delays  boundary  layer  reattachment  which  greatly  increases  the  size  of  the 
separated  region.  This  causes  the  shock  waves  to  move  forward  and  upwards  away  from  the  wall. 
The  resulting  shock  system  is  differently  shaped  compared  with  the  zero  pressure  gradients  case. 
It  follows  therefore  that  results  for  normal  shock  wave-boundary  layer  interactions  without  a 
post  shock  pressure  gradient  are  poor  guides  to  interactions  in  most  practical  flows  which  involve 
post  shock  pressure  gradients. 


1  Values  of  7>  are  unreliable  here  as  measurements  with  the  temperature  probe  facing  down¬ 
stream  into  the  reversed  flow  were  not  attempted. 
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A  plausible  explanation  of  boundary  layer  behaviour  in  flow*  with  normal  shock  waves 
can  be  made  in  terms  of:  the  shock  wave,  the  post  shock  adverse  pressure  gradient  and  vortices 
shed  from  the  interaction  region.  The  shock  wave  dominates  flow  development  in  its  immediate 
region  but  although  its  effects  (separation  and  a  rapid  rise  in  displacement  thickness)  are  drastic, 
they  are  also  shortlived.  Similarly  an  adverse  pressure  gradient  tends  to  increase  the  displace¬ 
ment  thickness  and  decrease  the  skin  friction  of  the  layer  but  these  effects  can  be  significantly 
modified  by  mixing  caused  by  longitudinal  vortices.  For  flows  with  a  shock  wave  the  longitudinal 
vortices  quickly  restabilize  the  boundary  layer  after  reattachment,  causing  the  skin  friction  to 
overshoot  its  undisturbed  value.  For  flows  with  a  shock  wave  and  a  subsequent  adverse  pressure 
gradient,  longitudinal  vortices  continue  to  move  the  layer  away  from  the  condition  of  a  corres¬ 
ponding  shock  free  reference  layer.  The  vortices  appear  to  be  very  long  lived,  as  this  process 
continues  to  distances  large  compared  with  distances  of  practical  interest. 

The  temperature  results  show  that  the  shape  of  the  mean  temperature-velocity  profiles 
depend  on  the  sign  of  the  pressure  gradient  and  the  local  heat  transfer  rate.  Profiles  can  be 
quickly  and  markedly  distorted  by  unusual  local  flow  conditions  such  as  a  nearby  separation 
or  reattachment  point  or  by  backflow  within  a  separation  region.  In  adverse  pressure  gradient 
flow,  the  temperature  velocity  profiles  lie  between  the  Crocco  formula  and  the  free  stream 
stagnation  temperature.  In  favourable  pressure  gradient  flow  they  lie  between  the  Crocco  and 
the  modified  Crocco  formulae. 
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FIG.  1  EXPERIMENTAL  RIG 
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FIG.  2  SCHLIEREN  SPARK  PHOTOGRAPH  OF  NORMAL  BIFURCATED  SHOCK  WAVE. 
The  wall  is  approximately  1/1  Oth  of  the  window  diameter  below  the  bottom 
of  the  window. 
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FIG.  4  SEPARATED  REGIONS  BEHIND  PRISMS  ATTACHED  TO  THE  WALL 
IN  A  TURBULENT  BOUNDARY  LAYER 
(Schofield  &  Logan  (1983)  ) 


FIG.  5  WALL  PRESSURE  DISTRIBUTIONS 


G.  6  MACH  NUMBER  DISTRIBUTIONS 


FIG.  7  WALL  PRESSURE  DISTRIBUTIONS  SHIFTED  TO  SHOW  SIMILARITY. 


FIG.  8  PRESSURE  GRADIENTS. 


FIG.  9  WALL  PRESSURE  DISTRIBUTIONS  NEAR  SHOCK  WAVE. 
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FIG.  11  STATIC  PROFILES- LAYERS 


MEAN  STREAMLINES  ON  CENTRELINE  IN  VICINITY  OF  THE  Note  vertical  scale  four  times  largerrfriaoihprizontal  scale 

INTERACTION.  Units  for  ^  are  kgm/m.s. 


FIG.  13  COMPARISON  OF  MEAN  STREAMLINES  FOR  FLOWS  WITH  Hatched  boundary  shows  extent  of  supersonic  torque. 

AND  WITHOUT  A  POST  SHOCK  PRESSURE  GRADIENT  Broken  line  indicates  edge  of  boundary  layer. 


FIG.  15  DISPLACEMENT  THICKNESS  GROWTH  WITH  AND  WITHOUT 
A  SHOCK  WAVE  -  HIGH  PRESSURE  GRADIENT  REGION. 


DISPLACEMENT  THICKNESS  GROWTH  IN  THE  MIXING  REGION 
WITH  AND  WITHOUT  A  SHOCK  WAVE. 
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FIG.  17  SHAPE  FACTOR  VARIATION  IN  THE  REGION  NEAR  THE  SHOCK  WAVE. 


FIG.  18  SHAPE  FACTOR  VARIATION 


FIG.  19  CLAUSE R'S  SHAPE  FACTOR 


FIG.  20  COLES’ 


FIG.  23  WALL  TEMPERATURE  DISTRIBUTION  -  LAYER  S 


OlOT) 


<u  <o 


RG.  24  ENTHALPY  THICKNESS  DISTRIBUTION  -  LAYER  S 


TEMPERATURE  DISTRIBUTION  -  LAYER  L 


WALL  TEMPERATURE  DISTRIBUTION  -  SUPERSONIC  FLOW  IN 
FAVOURABLE  PRESSURE  GRADIENT. 


1  0.2  0.3  0.A  0.5  0.6  0.7  0.8 


FIG.  27  TEMPERATURE  PROFILES  -  SUPERSONIC  FLOW  IN 
FAVOURABLE  PRESSURE  GRADIENT. 
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FIG.  27  (Concluded) 
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FIG.  29  TEMPERATURE  PROFILES  -  LAYER  S 
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FIG.  29  (Continued) 
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APPENDIX 

Mean  Velocity  Profiles— Plotted  on  Wall  Similarity  Co-ordinates,  u/Ui  versus  logioy  CA/v. 
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